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Abstract. The relation between fixed point and orbit count sequences is investigated from 
the point of view of linear mappings on the space of arithmetic functions. Spectral and 
asymptotic properties are derived and several quantities are explicitly given in terms of 
Gaussian binomial coefficients. 



1. Introduction and general setting 

Each map of an arbitrary set into itself gives rise to sequences a = (a n ) nS N and c = (c„) n eN 
that count the number of periodic points with period n and the orbits of length n, respectively. 
In many interesting cases, a n and c n are finite for all n £ N, turning a and c into sequences 
of non-negative integers. They are then related by Mobius inversion as 

(1) a n = J2 dc d and = - J^/i a d . 

d\n d\n 

Here, [i denotes the Mobius function from elementary number theory [lj. Puri and Ward 
|10j study these relations as paired transformations between certain non-negative integer 
sequences, with special emphasis on arithmetic and asymptotic aspects. It is the purpose 
of this short note, which is based on [8], to extend these transformations to the space of 
arithmetic functions, on which they act as linear operators, and to investigate their spectral 
properties and their asymptotic behaviour under iteration, thus generalising results from |10j . 

Let A = C N denote the set of arithmetic functions. Its elements can be thought of either 
as functions /: N — > C,n i— > /(n), or as sequences of complex numbers (f n )neN, hence the 
notations f(n) and f n are used in parallel; see pQ for details and general notation. Clearly, 
A is a complex vector space with respect to component- wise addition and the usual scalar 
multiplication. Additionally, it is equipped with the structure of a commutative algebra 
(A, +,*) with unit element, where * denotes Dirichlet convolution, defined by 

(/*<?) w =e <?g; 

d\n 

for all n £ N. The unit element with regard to this product is 5\ = (1, 0, 0, . . .). When / £ A 
is invertible, the notation / _1 refers to its Dirichlet inverse, defined by = / _1 *f = 8\. 

Furthermore, / • g denotes the pointwise product where (/• g){n) = f{n)g(n) for all n £ N. 
Motivated by the transformations (JTJ) , we define 'fix' and 'orb' as mappings from A into itself 
by a i— > a', with 

(2) a' n = fix{a) n = ^2da d and a' n = ovb(a) n = - £ //Q) a d , 

d\n d\n 
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respectively, for all n S N. Denoting by N the arithmetic function n \— > n, fix and orb can be 
written as 

(3) fix(/)=JV ■(!*/) and orb(/) = 1 • ( M * /) . 

This is particularly useful when considering the Dirichlet series of fix(/) and orb(/). In 
general, the Dirichlet series Df of / € A is defined by 

dm - 1 /(m) 



m" 

m=l 



which is viewed as a formal series at this stage. This generating function gives rise to the 
(formal) identities pQ 

D .f+9 = D f+ D 9 and D f*9 = D f D 9- 

Using these in connection with ([3|), a straightforward calculation shows 
Lemma 1.1. Let D a (s) denote the Dirichlet series of a £ A. Then, 

D &xia) (s) = ((s)D a (s-l) and D orHa) (s) = - * + £> tt (a + l). 

where £ is Riemann's zeta function. □ 

For a discussion of the identities from Lemma ll. II in a different context, and for interesting 
recent development on Dirichlet series as orbit counting generating functions, we refer the 
reader to [3J Sees. 2 and 10]. We now turn our attention to the structure of fix and orb as 
linear mappings on A. 



2. Spectral properties 

When studying fix on A, it is natural to ask for eigenvalues and eigenvectors. Let x be an 
arithmetic function that satisfies the fixed point condition Hx(x) n = Y^d\n dx d = x n for all 
n € N. Solving for x n yields the recursive relation 

(4) x n = j— dx d 

n>d\n 

for n > 1, so that all x n are uniquely determined from x 1 . The x n depend linearly on x\, 
which is the only degree of freedom in solving (|4|), wherefore the eigenspace for the eigenvalue 
1 is one-dimensional. Setting x\ = a\ = 1, we denote the resulting solution of @ by 
a = (a 1? o 2 , « 3 , . . .). Since fix maps multiplicative sequences to multiplicative sequences, it is 
an interesting question to what extent this is reflected by a. 

Lemma 2.1. Let a satisfy fix(a) = a, with a 1 = 1. The eigenvector a is multiplicative, so 
that a nm = a n ■ a m for all n, m > 1 with (m, n) = 1. 

Proof. In order to show multiplicativity by induction on the index, we first note that the 
claim is true for trivial decompositions such as n ■ 1, because a l = 1. In particular, it is then 
true for all natural numbers n < 5 = 2 ■ 3 — 1. Let (m, n) = 1 with mn > 1 and suppose the 
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ak satisfy the multiplicativity property for all k < ran — 1 (so that a k = a d a d , for all d, d! 
with (d, d') = 1 and dd' = k). This implies 



d|mn d=d m d n 
d<mn dm\m, d n \n 

d m d„<mn 



d|m d'|n 



This argument shows the induction step from ran — 1 to ran. In general, the step from N to 
N + 1 is trivial when N + 1 is prime, and of the above form otherwise, which completes the 
induction. □ 

For the remainder of this section, a always denotes the unique solution of fix(a) = a with 
ai = 1. The recursion in @ runs over all proper divisors of n and yields an explicit expression 
for a p r with p prime and r € N. By multiplicativity, this extends to a closed formula for a n . 

Proposition 2.2. Let n > 1 with prime decomposition n = p 1 ^ 1 ■ . . . -p r s s . Then, the n-th entry 
of the eigenvector a of Lemma 12.11 reads 

s r k 

a "=nn7— 7 

k=ie=i l p k 

Proof. Observing that a p = induction on the exponent r leads to a pr = Yie=i(^ ~P^)~ 1 > 
from which the statement follows by Lemma 12. 11 □ 

The sequence of denominators is now entry A153038 of the Online Encyclopedia of Integer 
Sequences (OEIS) [11]. In particular, \a n \ < 1 for all n, and Proposition 12.21 suggests that the 
a n become 'asymptotically small' in some sense, as n — > oo. Indeed, with 



i q = {feC m \J2\fn\ q <oo}, 



n=l 

one obtains the following result. 

Proposition 2.3. The eigenvector a is an element of i 1+£ for all e > 0. 

Proof. Substituting p k — 1 = (p — 1) Y^j=aP' i n the denominator of |a p r| leads to 

r r k—1 

jj(p* _ i) = (p _ iy yi y^p 1 > (p - i) r f ELi (fc_1) = (p- i)v (r - 1)/2 . 

k=l fc=l j=0 

Since (p — l) r p r ( r_1 )/ 2 > p r for r > 2 or p > 2, the inequality |ey | < holds in these cases. 
Besides, one has \a p \ = = • = for every prime p. The remaining term a 2 2 = 352 
introduces an additional factor of c = | in the estimate of \a n \, to account for the case when 
2 2 is the highest power of 2 that divides n. For n = p^ 1 ■ . . . ■ p r s s , this gives 

c tt Pi 



I On I < - n 



{i\n=i} 11 
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This product can be related to Euler's ^-function, 



p\n ' p\n p\n 

p 2 \n 



giving n {l |r l= l} pf-L — din)- Usm g 



lllXL Slip —— €^ 

n-,00 (f>(n) log(log n) 



where 7 = linin^oo (l + ^ + • • • + r — log re) denotes Euler's constant, compare [H Thm. 328], 
we obtain 

C • log(log(n)) 

On < — 



n 

for some constant C, which implies Y^m=i l a n| 1+£ < 00. □ 

In particular, a £ and (numerically) one finds X^Lil a n| 2 ~ 2.99635 < 3. Since J2 p | 
diverges and j a | = — > - , it is clear that the statement of Proposition 12.31 is not true for 

£ = 0. VP 

Lemma 2.4. Let D a (s) be the Dirichlet series of the eigenvector a, with s = a + it. 

(i) D a (s) converges absolutely in S = {s € C | a > 1}. 

(ii) linio-^oo D a {a + it) = 1, uniformly in t € R. 

Proof. We employ the methods of [TJ Chapter 11]. The first claim follows from the absolute 

convergence of ((a) for a > 1, because £) m < E m ^pr < C(»- 
For (ii), observe that a% = 1, wherefore one has 



00 

- 11 ^ E ^ = - 1, 



m=2 



which tends to (independently of t) as a — > 00. □ 



Using the fixed point relations fix(a) = a and orb(a) = a, the first equation of Lemma 1 1.1 
becomes a 'recursion' that can be solved for an explicit expression in terms of the Riemann 
zeta function. 

Theorem 2.5. The Dirichlet series of the fixed point a satisfies 

1=1 ^ K ' p e>i v 
where the infinite product converges on the set S from Lemma 12.41 

Proof. Using Lemma [LTj one inductively obtains D a (s) = TW^n^ f° r arbitrary m G N. 
For m — ► 00, the numerator converges to 1, since 



D a (s + m) = a 1 + Y: ~ <a 1 + ±D a { & 



n° n' 

n>2 



Due to the asymptotic behaviour Q(o + it) = 1 + 0(1/2°") for a — > 00, the product in the 
denominator of the resulting expression D a (s) = YYeLi £(s+e) conver g es as well, giving the 
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first part of the claim. The second equality follows from the Euler product representation 
£(s + i)^ 1 = rip(l — p~^ s+ ^) for a > 1, in which the order of taking products over p and I 
may be exchanged due to absolute convergence. □ 

Remark. Since a 1 7^ 0, a admits a Dirichlet inverse. Its terms can be calculated via to the 
general recursion formula, compare [TJ Thm. 2.8]. The Dirichlet inverse satisfies Dj,_ 1 = \/D^ 
wherefore b := a^ 1 has the Dirichlet series D b = Jlf>i C( s + ^)- Determining the coefficients 
of this series leads to an analogue of Proposition 12. 2 ( the n-th entry of b is 

• mi , 11 



=ifc=i 



f = \u'lM- 1 ) 



with n = p^ 1 ■ . . . ■ p r s s as before; see [8] for details. The first few terms of a and b read 

1,-1, ,-, ,-, , , • • -1 and b = (l, 1, — , . 

' 2 3 4 2 6 21 / V 2 3 4' 2 6 21 

Although it appears that a and b share the same denominator sequence, this is not true when 

the terms are represented as reduced fractions, because cancellation can occur for the b n . The 

first instance of this phenomenon is &12 = 3, whereas a\2 = ~\- 

Returning to the task of identifying all eigenvalues and their corresponding eigenspaces, 
we consider the equation x = A • fix(x) for A ^ 1 or, equivalently, 

(X-n)x n = A dx d 

n>d\n 

for all n S N. Since fix(x) 1 = x l by definition, the first term of such an eigenvector has to be 
0. In general, if x d = for all d\n with d < n, the n-th term x n will vanish as well, unless 
A = n. Thus, the first non-zero term is x m where m = A, and all eigenvalues are thus natural 
numbers. For fixed 1 < m E N, the recursion analogous to the case m = 1 given in (JH) is 

(5) x n = — - — V dx d (for n^m), 

m — n ^-^ 

n>d\n 

which shows that x n 7^ if and only if n = km for some k > 1. The recursion for the 
non-vanishing terms simplifies to 

v ' k>d\k k>d\k 

for k > 1, where the only free parameter is now the first non-zero term x m , indicating a one- 
dimensional eigenspace. Choosing x m = a 1 = 1, the resulting eigenvector a*-" 1 ) with eigenvalue 
m can be expressed in terms of a = a^- 1 ': 



(6) a^ m \n) 



a 



0, otherwise. 



In particular, G £ 1+£ for all e > and all m > 1, where a^" 1 ^ and a clearly have the same 
norm. Being eigenvectors for different eigenvalues, the are linearly independent. The 
next proposition states that they even generate the space of arithmetic functions in the sense 
of formal linear combinations. 
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Proposition 2.6. The eigenvectors of fix form a basis of A, so that \ m E N) c = A, 

and the representation of f E A as a formal linear combination of the is unique. 

Proof. Let / E A, a = and b := (a) -1 , and set a = f *b. Then, a*a = (f*b)*a = f. 
On the other hand, one has 

( a * a ) (k) = £a d a (d) (fc) = £ a m a (m) (fc) 
d|fc m>l 

for all k E N, so / E | m E N) c . Uniqueness of the linear combination follows from the 

uniqueness of the Dirichlet inverse. □ 

The spectrum of an operator T : A — > A is defined by 

spec(T) := {A E C [ T — A id is not invertible}. 

Note that fix and orb could be investigated as operators on I 2 , restricting their domains 
to appropriate invariant subspaces. But since we need a more general setting later when 
considering their behaviour under iteration, we work with the spectrum as a purely algebraic 
notion. 

Theorem 2.7. The spectra of fix and orb are 

spec(fix) = N and spec(orb) = |— JVgN 

Proof. Obviously, fix(/) = Xf implies / = if and only if A E C \ N, so fix — A id is injective 
for these A. Surjectivity can be seen by constructing a preimage / E A of an arbitrary g E A 
under fix — Aid. To this end, let /i = and, having already defined fx, . . . , f n -i for some 
n > 1, set 



1 



n>d\n 

When A ^ N, the resulting arithmetic function / is a term-wise well-defined element of A with 
fix(/) — Xf = g. In summary, fix — Aid is invertible if and only if A E C \ N. The spectrum 
of orb can be derived from that of fix by the equivalence 

fix(aW) = Na W ^ orb(oW) = la^ 

which follows from fix o orb = orb o fix = id. □ 

3. Iterations of the operator fix 



Puri and Ward [10] raised the question of the nature of the orbits under fix. In particular, 
the asymptotic properties of such orbits are sought. In our terminology, this is to ask how 
elements of A behave under iteration of the operator fix. Starting from 5\ = (1, 0, 0, . . .), the 
first iterates are sequences A000007, A000012, A000203, A001001 and A038991 - A038999 of 
the OEIS, see [UllCLl]. Here, we start from the sequences 6 m , m > 1, defined via 5 m (n) = 8 m ^ n 
for all n, and give term-wise exact representations of fix n (5 m ), from which the corresponding 
quantities of general starting sequences can later be constructed. 

Let y(°) = Si and := fix(y( n - 1 )) for n > 1, so that 

yW(m) = (fix(l/ (n_1) )) m = (fix n (2/ (0) )) m 

holds for all m > 1. 
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Proposition 3.1. Let p be a prime and r G N. Then 

r-l 

y (n) (p r ) 



n + r — 1 
r 



i „n+r—\—i 

-nV 

P i=0 



p 



i+1 



where ["] denotes the Gaussian or q -binomial coefficient. 

Proof. The claim can be verified by induction in n. By definition of fix, we have 

r 



y 



d\p r 



k=0 



(7) =1+ P lt 

P k=2 

Applying Pascal's identity for Gaussian binomials 



m 




m — 1 




m — 1 


t 


r 

= q 


t 


+ 


£-1 




1 




q 





to the fc-th summand in ([7]) yields 



n + k — 1 




n + k 




n + k — 1 


k 


p 


k 


p 


fc-1 



Plugging this and the relation [™ ] = ^2^=0 P k ^ n *° ©> we obtain 



y (n+1 \p r ) 



n+l 
1 



P k=2 



n + k 
k 



n + k — 1 
fc - 1 



n + r 
r 



which completes the proof. 

Corollary 3.2. For 1 ^ m € N itfii/t prime decomposition m = p^p^ 2 ■ ■ ■ ■ ■ p r k k , one has 



□ 



i=i i=0 



Proof. It is easy to check that fix, and hence all iterates fix n with n > 1, preserve multiplica- 
tivity. Thus, in view of Proposition 13.11 the claim follows from 6i being multiplicative. □ 

Remark. A different approach to calculating y^ exploits the fact that its Dirichlet series 
D ^(s) is known from the context of counting sublattices of Z n . Let 

h n {m) = \{A I A is a sublattice of Z n of index m}| 

and Dh n (s) = Y^ m >l m» • Proposition A.l of [2] states that 



h n (m) 



£ and D hn (s) = {(s)-D hn _ 1 (s-l) 



df ...■d n =m 



where the first sum runs over all n-tuples (d\, . . . , d n ) of positive integers with di-...-d n = m. 
In particular, D hi (s) = ((s) = D ^(s). On the other hand, Lemma fTTI gives D {n+1) (s) = 
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C(s)-D (n) (s — 1), so -D („) (s) = D hn (s) for all n > 1 and thus h n (m) = y^ n '(m) for all n, m > 1. 
For m = p^ 1 • . . . • pT* , the coefficient h n (m) can be written as 



k r j n+i-l _ -1 k 



n + rj — 1 



p.; 



as shown in [5] and [12] , leading to the result of Corollary 13,21 This also provides a concrete 
interpretation of the iteration sequences mentioned earlier. 

Similarly to calculating fix-eigenvectors for m > 1, the sequences arising from fix- 
iterations on starting sequences 5 m with m > 1 can be related to the case m = 1. From 
(fix(<5 m )) fc = m • 4,jm fo r some j € N, one can conclude inductively 



(fix n (<W), 



^•y (n) (£) 

0, otherwise. 



Since arbitrary elements from A can be written as complex linear combinations of the <5 m , we 
find the following behaviour of general arithmetic sequences under fix- iterations. 

Lemma 3.3. Let f = (fi, f 2 , ■ ■ ■) G A. Then, for all M > 1, 

M 

( fi * n (/))M = £a (fi^(^)) M = £/ fc *V n) (f )• 

fc=l fc|Af 

j4s a convolution product, this is 

fix" (/) = (/ -JV n ) /oralliiGN, 
where N n denotes the arithmetic function defined by m i— > m n . 

Proof. Note first that, for all M > 1, the value fix(/)j\/ depends on the / n with n|M only. Each 
single term of fix(/)^ can thus be obtained by applying fix to the finite linear combination 
Efc=i fkb~k, resulting in a component- wise well-defined arithmetic function. The last identity 
follows from the definition of the convolution product. □ 

Given the unbounded spectrum of fix, it is not surprising that the image sequences fix n (/) 
of a non- negative / component- wise tend to infinity when n — ► oo. Nevertheless, the 'quotient 

fix n (/)M ) nave a well-defined convergence behaviour for each M as follows. 

Theorem 3.4. Let f be a non-negative arithmetic function and let M 6 N be such that 
(fix r (/)) M ^ for some r 6 N. Then, 

(fix" +1 (/)) M 



iim \_ = M. 



r™ (fix»(/)) 



Proof. Expanding the quotient in question yields 

( fi * n+1 (/)) M _ ^ d (fix"(/)) d _ M+ _ d (fi*»(/)) t 



( fiX?l (/))M (^(/))m M^M ( fiXn (^)M 
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The number of terms in the last sum being independent of n, it suffices to show that, for each 
divisor d of M , the corresponding summand approaches as n — ► oo. Let d 1 = 1, . . . , d R = d 
denote the divisors of d. Then, using Lemma 13.31 one gets 



(fix«(/)) M ^ Ek\ M fkk n yW{%) ~ U d? UyWiM/di) jr{ yW(MM)' 

This upper bound can be seen to converge to by splitting numerator and denominator of 
each summand in the last sum into the factors that arise from the prime decomposition of 
d/di and M/di. For r > s, the quotient 

y {n) {pS) ^ s-1 1 _ p n+s-l- l r-l 1 _ pi+1 
^ ' y( n ^(p r ) 1 — p n + r ~ 1 ~ i J»A 1 — pn+r—l—i 

converges to as n — > oo, since the first product tends to p s ~ r > 0, while the second tends 
to 0. Being composed of factors such as ©, the quotient (fix n (/)) d /(fix ri (/)) M tends to 
for all d\M as n — > oo, which proves the claim. □ 

Remark. Unlike fix, the operator orb's spectrum is bounded by from below and by 1 from 
above. One would expect orb n (/) to converge to an element of the eigenspace for the largest 
eigenvalue (A = 1), hence to a scalar multiple of a^ 1 ' = a. Indeed, employing the eigenvector 
basis representation / = Yl m >i a m a ■> one obtains 

M M 
lim (orb n /) M = lim V a m orb n (a (m) )Af = lim a m n a M = a l a M 

n^oo n^oo ' — » n~ >oo £ — • m 

m=l m=l 

for all M e N. In other words, orb n (/) ra ^°° > f x a, for all / G A 

4. Concluding remarks 

The eigenvector a has an interesting consequence in the setting of dynamical (or Artin- 
Mazur) zeta functions, where one considers the generating function 

(OO j, \ oo 

E — z ™ ) = II( 1 - zm ) _orb(/)m ' 
m=l / m=l 

usually for arithmetic functions / that count (isolated) fixed points of a discrete dynamical 
system, see [3] for background material. Using a together with a = orb(a), one finds the 
product identity 



\m=l / m=l 

or, equivalently, the series identity 

oo oo 

Z — J nm — ^ 



m 

m=l m=l 



with absolute convergence for \z\ < 1. It would be interesting to know whether this gives rise 
to another interpretation of the meaning of the eigenvector a in this setting. 
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